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Abstract. Cosmology provides an excellent laboratory for testing various aspects of neutrino physics. 
Here, I review the current status of cosmological searches for neutrino mass, as well as other properties of 
>«-2 , neutrinos. Future cosmological probes of neutrino properties are also discussed in detail. 
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1. Introduction 

The last few years has seen a dramatic increase in the amount and precision of cosmological data. Very 
precise measurements of the cosmic microwave background has been performed by the WMAP satellite [1] 
and a number of different ground based experiments, most notably the ACBAR [5] and QUAD [3 telescopes. 
The large scale distribution of galaxies has been measured by the Sloan Digital Sky Survey (SDSS) project 
which released its final data in 2009. A number of other experiments have probed the formation of cosmic 
structure using other techniques, some of which are still relatively new. Altogether, the cosmological data 
is now precise enough that it is possible to probe various aspects of particle physics using cosmology. For 
example the density of dark matter has been measured quite precisely and since the dark matter particle is 
likely to be associated with physics at the TeV scale, i.e. with the energy scale probed at the LHC, cosmology 
provides an important second window to physics at this energy scale. 

However, perhaps the best example of the interplay between cosmology and particle physics is in neutrino 
physics. Here, precision cosmology can be used to probe questions traditionally investigated in laboratory 
experiments. While cosmology is at present not sensitive to the neutrino mass differences, it is highly 
sensitive to the absolute neutrino mass. Furthermore, this parameter is notoriously difficult to probe in 
laboratory experiments. During the past 7-8 years neutrino cosmology has been a source of rapidly growing 
interest and a large number of papers have been written on the subject. There are also several good review 
articles available in the subject (e.g. [HE]) focusing on various aspects of neutrino cosmology. In this review 
I focus mainly on the effects of neutrinos in cosmic structure formation, with particular emphasis on the 
observational probes which can be used to measure neutrino masses. 

Section 2 discusses the physics of neutrino decoupling in the early universe, including the effects of 
neutrinos on Big Bang nucleosynthesis. Section 3 briefly discusses the fact that during cosmic structure 
formation, long after neutrino decoupling, cosmic neutrinos are pure mass eigenstates and retain no flavour 
memory. In section 4 I review cosmic structure formation in some detail including both linear perturbation 
theory as well as neutrinos in non-linear structure formation. Section 5 is a review of the current constraints 
on neutrino mass, and section 6 contains a discussion of the current bound on the total energy density in 
neutrinos or other light, weakly interacting species. Section 7 discusses various future observational probes 
as well as their potential and possible systematics, and section 8 contains a brief discussion of the cosmic 
neutrino background anisotropy. Finally, section 9 contains a brief discussion and conclusion. 

1.1. Neutrino mixing 

Neutrino oscillation experiments have shown beyond any reasonable doubt that neutrinos have non-zero 
masses. Furthermore, two distinct mass differences have been measured in different types of experiments, 
perfectly compatible with the assumption that there are three distinct mass eigenstates, vm, corresponding 
two the three known flavour eigenstates, vp. These states are connected via 

v F = Uv M , (f) 

where U is a 3 x 3 unitary matrix. For Dirac neutrinos the mixing matrix contains 4 free parameters, the 
three rotation angles, and one complex phase, S. For Majorana neutrinos there are two additional phases 
which cannot be rotated away. These Majorana phases have no influence on oscillation observables, but are 
important for example in neutrinoless double beta decay (see e.g. [SI [7] for details). 

The matrix U can be parameterized in a number of ways, but the most popular by far is in terms of 
the mixing angles and the phases. In the relativistic limit neutrino oscillation probabilities only depend on 
the masses of the three states via Am 2 , i.e. in addition to the elements of the mixing matrix oscillation 
experiments are sensitive to two Am 2 values. 

By now constraints on U and Am 2 come from a variety of different experiments, the most recent overview 
is given in Ref. [5] and provides the following constraints (at 3cr) 

0i 2 =(34 ±3)° (2) 

02 3 =(43^)° (3) 

6 13 = < 12.5° (4) 
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5m\ 2 = 7.6±0.6x l(T 5 eV 2 (5) 

8m% x = ± (2.4 ± 0.4) x 10~ 3 eV 2 (6) 

(5cp = No constraint (7) 

As will be discussed later, apart from light element abundances, cosmological observables in most cases 
depend little on flavour and conversely cosmology provides little information on the mixing structure of 
neutrinos. 

1.2. Neutrino masses 

While oscillation experiments are by far the most useful probe of the mixing angles and the mass differences, 
Am 2 , they are hardly useful for probing the absolute neutrino masses, m,. To measure these parameters 
there are currently three useful paths. One possibility is to look for neutrinoless double beta decay which 
can occur because the emitted left handed neutrino contains a negative helicity component with relative 
amplitude mjE, which can then be absorbed as an antineutrino, provided that neutrinos are Majorana 
particles. The amplitude for this process is proportional to |^- U^mj\ (see e.g. [7] for for details). This 
particular form comes from the fact that each neutrino vertex is proportional to U e j and only the "opposite" 
helicity component contributes. Note that the sum retains the phase structure so that phase cancelation is 
in principle possible. 

The actual parameter measured in any neutrinoless double beta decay experiment is the half-life T\i% 
for some isotope which is related to mpp — | ^ . U^rrijl via the relation 

J_ =G 0,| M 0,|2 TO 2 (g) 

1 1/2 

where G " is a phase-space factor and |M "| the nuclear matrix element squared. Translating any 
measurement of a finite lifetime to an effective neutrino mass therefore involves an uncertainty related 
to the calculation of the nuclear matrix element (see e.g. [T21 [T3] for a discussion of this point). The best 
upper bound on mpp from double beta decay currently comes from the Heidelberg-Moscow (HM) experiment 
and is is mpp < 0.27 eV (90% C.L.) [§1[T3]. 

However, there is also a claim of a positive signal for the decay in a different analysis of the (HM) 
experiment [TS1 HH1 HZ]. The upcoming GERDA [TU] and EXO [H] experiments which are currently in 
commissioning will have the capability to either confirm or rule out this claim. 

Another option is to look for the kinematical effect of a non-zero neutrino mass in ordinary beta decay. 
This is challenging, not because the process is forbidden for zero neutrino mass, but because the neutrino 
mass only produces a noticeable effect very close to the endpoint of the electron spectrum where there are 
very few events. In the beta decay vertex, an electron neutrino is emitted. However, since the subsequent 
measurement involves the electron energy it is more convenient to view the beta decay as three separate 
possible processes which contain either v\ , v 2 or v^ . The beta decay rate is the incoherent sum of these three 
processes, and the rate therefore involves ^2 U^Uej, i.e. it does not retain any Majorana phase information. 
The shape factor for the electron decay in the presence of more than one massive neutrino is given by 

S{E e ) ex (Q - E e ) J2 U *ej U ej ^{Q-E e ) 2 -m 2 . (9) 



As long as the energy resolution of the experiment is significantly worse than the mass splittings this can be 
written as 



S(E e ) oc (Q - E e ) (Q - E e ) 2 - J2 U; 3 U ej m 2 , (10) 

i.e. the spectrum distortion can be described by a single effective mass. The current upper bound on this 
effective mass is 

1/2 



mp 



Y; U *ejUe 3 m 2 A <2.3eV (11) 
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at 95% C.L. from the final analysis of the Mainz experiment [T3]. Starting in early 2012 the KATRIN 
experiment jT8l [19] will improve this sensitivity by an order of magnitude to 0.2 eV. 

The final method is to look for the kinematic effect of neutrino masses in cosmological structure 
formation. As will be discussed in detail in the next sections, a non-zero neutrino mass means that 
neutrinos contribute to the present density of dark matter. However, since their thermal history is very 
different from that of cold dark matter they have a very distinct signature on cosmic structure. As a first 
approximation, cosmology is sensitive to the total energy density in neutrinos, which for non-relativistic 
neutrinos is simply proportional to the sum of all neutrino mass eigenstates, X) m j- However, this is 
true if the sensitivity of cosmological data to the neutrino mass is poor compared with the internal mass 
splittings, i.e. if a(m) > ^ m v . Some proposed cosmological structure surveys in the coming decade will be 
sensitive to masses as low as 0.03-0.05 eV, which means that the effect of massive neutrinos cannot simply 
be approximated as a sum of the involved mass eigenstates. This possibility will be discussed later. 

The complementarity of these three ways to probe neutrino masses has been discussed in some detail 
in the literature. While they all probe the neutrino masses, they do so via different combinations of the 
mass states with the mixing matrix elements. They are also prone to completely different systematics, and, 
just as important, to different modeling uncertainties. For example, the presence of right-handed currents 
will affect both the neutrinoless double beta decay rate and the shape of the beta decay spectrum, but have 
no impact on the cosmological mass measurements. It is therefore entirely possible that an experiment like 
KATRIN measures a spectrum distortion, while cosmology provides no signature. 

2. Neutrino decoupling 

Before going on to discuss the impact of neutrinos on structure formation I will briefly review the standard 
picture of neutrino decoupling in the early universe. 

In the standard model neutrinos interact via weak interactions with charged leptons, keeping them in 
equilibrium with the electromagnetic plasma at high temperatures. Below T ~ 30 — 40 MeV e + and e~ are 
the only relevant particles, greatly reducing the number of possible reactions which must be considered. In 
the absence of oscillations neutrino decoupling can be followed via the Boltzmann equation for the single 
particle distribution function [20] 

m~ Hp dp- = Cco11 ' (12) 

where C co n represents all elastic and inelastic interactions. In the standard model all these interactions are 
2 •(-> 2 interactions in which case the collision integral for process i can be written 

d 3 p 2 d 3 p 3 d 3 p 4 



a 



coll 



iifl) -2E~J ' 



2E 2 (2tt) 3 2E 3 (2tt) 3 2£ 4 (2tt) 3 
x (2ir) 4 6 4 ( Pl +P2-P3+P4)A(fij2j 3 ,h)S\M\l 2 _ 3 ^, (13) 

where S\M\\ 2 ^ 3A i is the spin-summed and averaged matrix element including the symmetry factor S = 1/2 
if there are identical particles in initial or final states. The phase-space factor is A(/i,/2, /3, /*) = 

UU0- - /i)(i - h) - /i/ 2 (i - /s)(i - h). 

The matrix elements for all relevant processes can for instance be found in Ref. [21] (however, see also 
[5]). If Maxwell-Boltzmann statistics is used for all particles, and neutrinos are assumed to be in complete 
scattering equilbrium so that they can be represented by a single temperature, then the collision integral 
can be integrated to yield the average annihilation rate for a neutrino 

T= l -^{gl+gl)T\ (14) 



where 



sin 6w + (r + sin 9w) f° r v & 

9l+9R=\ f \ ■ (15) 

sin 4 0w + (— r + sin 2 9w) 2 for v^ r 
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This rate can then be compared with the Hubble expansion rate 

H = IMgl' 2 ^— (16) 

Mpi 

to find the decoupling temperature from the criterion H = T| T=T . From this one finds that Tr){v e ) — 2.4 
MeV, Toiy^r) — 3-7 MeV, when g* = 10.75, as is the case in the standard model. 

This apparently straightforward conclusion is complicated by neutrino oscillations. As has been shown 
in a number of papers the mass differences and mixing angles of the neutrino sector are such that all species 
are effectively almost equilibrated prior to neutrino decoupling. To a good approximation neutrinos can 
therefore be treated as a single species with an averaged coupling strength, decoupling slightly prior to the 
estimated decoupling temperature of electron neutrinos. 

A coupling temperature of around 2.5-3 MeV means that neutrinos decouple at a temperature which is 
significantly higher than the electron mass. When e + e~ annihilation occurs around T ~ m e /3, the neutrino 
temperature is unaffected whereas the photon temperature is heated by a factor (11/4) 1 / 3 . The relation 
T v /T~ f = (4/11) 1 / 3 ~ 0.71 holds to a precision of roughly one percent. The main correction comes from a 
slight heating of neutrinos by e + e~ annihilation, as well as finite temperature QED effects on the photon 

propagator m m OH H3 \M mi EH M M M ED GH 133 GH 125] • 

The extra energy deposited in neutrinos is usually defined in terms of an extra number of neutrino 
species, AN V , defined as 

AN U = ^-, (17) 

Pv,0 

where p v> o is the energy density in a single neutrino species assuming complete decoupling. 

The most precise calculation to date [35] has estimated that AA„ = 0.046, so the standard model 
prediction for the total energy density in neutrinos is N v — 3.046. As will be discussed later the difference 
of 0.046 is probably too small to be detectable even with future observational data. 

It should be noted here that N v is customarily used in cosmology to parameterize any additional 
relativistic energy density, not just neutrinos. Therefore N v is one of the cosmological parameters normally 
fitted in cosmological parameter estimation and any value significantly different from 3 could indicate the 
presence of new physics beyond the standard model. 

2.1. Big Bang nucleosynthesis and the number of neutrino species 

Shortly after neutrino decoupling the weak interactions which keep neutrons and protons in statistical 
equilibrium freeze out. Again the criterion H — T\ T=T can be applied to find that Tf reeze ~ 0.5.9, 

MeV ng. 

Eventually, at a temperature of roughly 0.2 MeV deuterium starts to form, and very quickly all free 
neutrons are processed into 4 He. The final helium abundance is therefore roughly given by 

1 + n n /n 
n n /n p is determined by its value at freeze out, roughly by the condition that «■«/%> |T=T froozc ~ 

g- {ra n —m p )/Tf„ M „ 

Since the freeze-out temperature is determined by g* this in turn means that g* can be inferred from a 
measurement of the helium abundance. However, since Yp is a function of both flf,h 2 and g* it is necessary 
to use other measurements to constrain ft^h 2 in order to find a bound on g*. Historically, this has been done 
by using the measured abundance of deuterium as a probe of the cosmic baryon density. However, at present 
the most precise determination of the baryon density by far comes from the cosmic microwave background 
observations. 

When the baryon density is fixed, Yp can be used to constrain g*. As will be seen later, observations 
of the cosmic microwave background also provide a constraint on g* which is quite stringent and the two 
types of observations can be used as a consistency check for the standard radiation dominated expansion 



Y P 



(18) 

T~0.2 MeV 
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in the region eV < T < MeV. Any discrepancy between the two values of <?* could in principle indicate 
non-standard physics such as the late decay of a massive particle. 

Usually such bounds are expressed in terms of the equivalent number of neutrino species, N v = p/p„ , 
instead of <7*. The exact value of the BBN bound on N v is somewhat uncertain because of systematic 
uncertainties involved in the Yp determination. 

Very recent analyses are those found in [36, 37 which indicate Yp = 0.2565 ± 0.001 (stat) ± 0.005 (syst) 
[3"6] and Yp = 0.2561±0.011 [37] respectively (see also [38.). The central values are consistent with N v ~ 3.7, 
i.e. a value somewhat higher than the predicted N v = 3.04. Even though this might point to non-standard 
physics, the possible systematics seem too large for any firm conclusions. However, it is intriguing that CMB 
observations currently point in the same direction, i.e. N v > 3. 

Another interesting parameter which can be constrained by the same argument is the neutrino chemical 
potential, £„ = [i v jT [39, %01l41[ H2]. At first sight this looks like it is completely equivalent to constraining 
N v . However, this is not true because a chemical potential for electron neutrinos directly influences the n— p 
conversion rate. Furthermore, it is crucial to take neutrino flavour oscillations into account when calculating 
bounds on the neutrino chemical potential. This point is discussed in more detail below. 

2.2. The effect of oscillations 

In the previous section the one-particle distribution function, /, was used to describe neutrino evolution. 
However, for neutrinos the mass eigenstates are not equivalent to the flavour eigenstates because neutrinos 
are mixed. Therefore the evolution of the neutrino ensemble is not in general described by the three scalar 
functions, /,, but rather by the evolution of the neutrino density matrix, p = ijnp' , the diagonal elements of 
which correspond to /,. 

Using the density matrix, p, the Boltzmann equation is replaced by (see e.g. |138p 



dp _ R dp_ 
dt dp 



M 8V2G F p 
2p 3m^ 



■C[p\. (19) 



The first term accounts for oscillations and the last term is the 0(G F ) collision operator. M. is the mass 
matrix in flavour space, i.e. M. = U'MU, with M = diag(m^,?7i|,m§). 

Eq. [in] was solved approximately for the case of standard neutrinos in Ref. [31] , and exactly in [35] . 

Without oscillations it is possible to compensate a very large chemical potential for muon and/or tau 
neutrinos with a small, negative electron neutrino chemical potential [39] . However, since neutrinos are 
almost maximally mixed a chemical potential in one flavour can be shared with other flavours, and the end 
result is that during BBN all three flavours have almost equal chemical potential [55] H3J 2H 25J 25] ■ This in 
turn means that the bound on i> e applies to all species. The most recent bound on the neutrino asymmetry 
from BBN comes from [69 

- 0.04 < & = 1M < 0.07 (20) 

for i = e, (J,, t|||. 

The bound assumes complete flavour equilibration during BBN, which with the measured mixing angles 
and mass differences is a fairly good approximation. It should, however, be noted that with some fine tuning 
of the initial conditions it is possible to have large lepton asymmetries while still producing the correct light 
element abundances [70] . Another possibility is that additional majoron type interactions may allow for 
large asymmetries to be present [7T] . 

In models where sterile neutrinos are present even more remarkable oscillation phenomena can occur. 
However, I do not discuss this possibility further, and instead refer to the review [S]. 

| It should be noted that this bound is strictly speaking only valid assuming that there is no extra relativistic energy density 
in other species present. 
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3. What is a cosmic background neutrino? 

An interesting question which is not often addressed in the literature is whether cosmic background neutrinos 
are mass or flavour eigenstates, or neither. Even if neutrinos are created as flavour states in charged current 
interactions, they stop interacting at the MeV scale and after this point they propagate freely. The neutrino 
created in a charged current interaction can be thought of as a single wavepacket consisting of a superposition 
of three mass states. These three wavepackets travel at slightly different velocity and at some point the 
original wavepackets will no longer overlap, i.e. the original flavour state has decohered and separated into 
distinct wavepackets related to mass states. 

The timescale for this decoherence phenomenon to occur can be estimated from the following argument 
(see [7J H7] for a more detailed discussion) . 

At creation, the neutrino wave packet should have a size of roughly 8 X — cAi co n, where Ai co n is the 
mean time between collisions for particles in the medium. To be conservative we assume the particles in 
the medium to interact only via the weak interactions which create the neutrinos. In practise electrons for 
example also have electromagnetic interactions which are much faster. Assuming only weak interactions is 
therefore an upper bound on the time between collisions and therefore also an upper bound on the size of 
the emitted neutrino wavepackets. 

At decoupling, by definition, the timescale for collisions is comparable to the Hubble time, At co ii ~ 
H^ 1 ~ mpi/Tp. At some later time the size of the wavepacket will have increased to 5 X = mp\/(TrjT) 
because of cosmic expansion. Now assume that the original wavepacket is a superposition of two 
wavepackets related to different masses, m\ and mi- The velocity difference between the two packets is 
Av ~ (m 2 , — m\)/p 2 = dm 2 /p 2 ~ 8m 2 /T 2 , and once the neutrino has traveled a distance of L ~ 5 X /Av 
the two wavepackets no longer overlap and the original state has decohered. The distance covered by the 
neutrino at some given time is the Hubble scale at that time, mpi/T 2 , again assuming that neutrinos 
are ultrarelativistic. This distance is longer than the distance needed for decoherence provided that the 
decoherence condition 

5m 2 > T 3 /T D , (21) 

is fulfilled. Note that this assumes radiation domination, in a matter dominated universe the condition is 
slightly different. At matter-radiation equality this corresponds to Sm 2 > 10~ 6 eV 2 which is fulfilled by 
both dm 2 1 and Sm^. We again note that this is a very conservative estimate, in practise charged leptons 
have much shorter collision timescales, and Sp is therefore correspondingly much larger than the pure weak 
interaction estimate. 

Thus, cosmic background neutrinos can be treated as exact mass eigenstates during the entire history 
of structure formation in the universe. 

This conclusion also leads to the question: What is the number density of neutrino states of mass state 
i in the universe? At high temperature long before neutrino decoupling, the neutrino density matrix is 
diagonal in flavour space, with the diagonal elements characterised by temperature and chemical potential, 
T and £ a . After decoupling the evolution becomes simple in mass basis, i.e. 

Uptf =p M = (22) 

in the absence of decoherence. Decoherence can be approximated by a damping term in Eq. 1221 such that 

p M = -K(p M - diag(p M )), (23) 

where if is a temperature dependent decoherence rate. In the absence of flavour asymmetry the evolution 
is unimportant, but in the presence of an asymmetry the final density of states 1,2 and 3 is given by 

rii = PM,ii(t = *dec) (24) 

with all off-diagonal components in the mass basis damped away. 
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4. Neutrinos in cosmic structure formation 

Neutrinos are a source of dark matter in the present day universe simply because they contribute to VL m . 
The present temperature of massless standard model neutrinos is T„,o = 1.95 A' = 1.7 x 10~ 4 eV, and any 
neutrino with m 3> T v $ behaves like a standard non-relativistic dark matter particle. 

The present contribution to the matter density of N u neutrino species with standard weak interactions 
is given by 

*-*" = N -i^hv < 25 > 

Just from demanding that £l„ < 1 one finds the bound [48J 09] 

/ 46 eV 
m v < -^— (26) 

More realistically, one could assume that all dark matter is in the form of neutrinos. Very loosely, the total 
matter density is fl m ~ 0.3, with approximately 0.05 in the form of baryons. A more realistic upper bound 
would then be £l„h 2 < 0.12, leading to 

a number comparable to the present upper limit from beta decay experiments. However, as will be seen 
below, observations of cosmic structure allow for a much tighter constraint on the neutrino mass. 

4-.1. The Tremaine-Gunn bound 

If neutrinos are the main source of dark matter, then they must also make up most of the galactic dark 
matter. However, neutrinos can only cluster in galaxies via energy loss due to gravitational relaxation since 
they do not suffer inelastic collisions. In distribution function language this corresponds to phase mixing of 
the distribution function |50| . Initially, the microscopic neutrino distribution function is given by the simple 
non-degenerate Fermi-Dirac distribution 

m = i^TT' (28) 

where T — To/a, i.e. the distribution function is preserved as an ultra-relativistic distribution function after 
neutrino decoupling in the early universe. 

By using the theorem that the phase-mixed or coarse grained distribution function must explicitly take 
values smaller than the maximum of the original distribution function one arrives at the condition 

/CG < fv,mssx = X (29) 

Because of this upper bound it is impossible to squeeze neutrino dark matter beyond a certain limit 50 . 
For the Milky Way this means that the neutrino mass must be larger than roughly 25 eV if neutrinos make 
up the dark matter. For irregular dwarf galaxies this limit increases to 100-300 eV [5TJ[52], and means that 
standard model neutrinos cannot make up a dominant fraction of the dark matter. This bound is generally 
known as the Tremaine-Gunn bound. 

Note that this phase space argument is a purely classical argument, it is not related to the Pauli blocking 
principle for fermions (although, by using the Pauli principle /„ < 1 one would arrive at a similar, but slightly 
weaker limit for neutrinos). In fact the Tremaine-Gunn bound works even for bosons if applied in a statistical 
sense [51j . because even though there is no upper bound on the fine grained distribution function, only a 
very small number of particles reside at low momenta (unless there is a condensate) . Therefore, although the 
exact value of the limit is model dependent, limit applies to any species that was once in thermal equilibrium. 
A notable counterexample is non-thermal axion dark matter which is produced directly into a condensate. 
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It should also be noted that in the original work by Tremaine and Gunn the coarse grained distribution 
was assumed to be Maxwell-Boltzmann. However, this does not seriously influence the bound. Kull et al. 
[72] derived a more general bound 

m»(^y /4 , (30) 

in which no specific assumption about the final distribution function is involved. 

A very interesting direct example of the Tremaine-Gunn bound was studied in [73] ■ Here, neutrino 
clustering in evolving cold dark matter halos was studied using the Boltzmann equation. The problem was 
made tractable by assuming no backreaction, i.e. that the CDM halos arc not affected by neutrinos. The 
results clearly show that there is a maximum in the coarse grained distribution function of 1/2, as expected. 
The calculation was extended to bosons in [53], where no such upper bound was found. In fact bosons can 
have an average density several times higher than fermions with equal mass inside dark matter halos, an 
effect which in principle might be used to distinguish fermionic and bosonic hot dark matter. 

Finally, the effect was studied in detailed iV-body simulations in [57] . This point will be discussed in 
more detail below. 

J^.2. Neutrino hot dark matter 

A much stronger upper bound on the neutrino mass than the one in Eq. (I26p can be derived by noticing that 
the thermal history of neutrinos is very different from that of a WIMP because the neutrino only becomes 
non-relativistic very late. 

The Boltzmann equation can generically be written as 

L[.f} = ^ = C[f}, (31) 

where L[f] is the Liouville operator. The collision operator C[f] on the right-hand side describes any possible 
collisional interactions. For neutrinos C[f] = after neutrino decoupling at T ~ 2 — 3 MeV. 
We then write the distribution function as 

f(x i ,q,n j ,v)=fo(q)ll + V(x i , q,n J3 v)], (32) 

where fo(q) is the unperturbed distribution function. For a fermion which decouples while relativistic, this 
distribution function is 

/o(<z) = [cxp(g/T ) + l]- 1 , (33) 

where To is the present-day temperature of the species. 

In conformal Newtonian (longitudinal) gauge the Boltzmann equation for neutrinos can be written as 
an evolution equation for "J/ in &;-space [54j 



T L ifl = it + lk ~^ + -n — 

jo or e dmq 



<p — ik-fj,ip 

q 



= 0, (34) 



where \x = n 3 kj. ip and <j) are the metric perturbations, defined from the perturbed space-time metric in the 
conformal Newtonian gauge [54] 

ds 2 = a 2 {r/)[-{l + 2i))dr? + (1 - 2<j))5 ij dx l dx j ]. (35) 

The perturbation to the distribution function can be expanded as follows 

oo 

* = EH) / (2Z + l)^ J P/(^)- (36) 
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One can then write the collisionless Boltzmann equation as a moment hierarchy for the ^('s by performing 
the angular integration of L[f] 

;dlnfo 



e 



*i=fcj(* 



*/ 



1 



din q 

. ek.dha.fo 
3q ding 

-(W?,_i-(J + l)ttl + i) , Z>2. 



(37) 

(38) 

'(2J + l)e v "~' -i v ~ ' "'-'^ > •--- ( 39 ) 

As long as neutrinos are relativistic the solution to this system of equations is damped for all modes 
inside the horizon. Examples of the exact numerical solutions, calculated using CAMB |55) are presented in 
Fig. [TJ When the neutrino becomes non-relativistic the gravitational source term (the k 2 ip term) becomes 
important and begins to feed the \l/'s, starting with the lowest multipoles. The higher multipoles are affected 
slightly later. Another important point is that \& starts to grow later for the smaller k value. This point will 
be discussed below. 




Figure 1. *I/;'s for 3 neutrino masses with momentum q/To = 3 as a function of the scale factor. The upper 
three panels are for k = 0.01/iMpc - 1 and the lower three panels for k = 0.1/iMpc . 
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By integrating the neutrino perturbation over momentum 

Kl Jdqq 2 ef (q) ' ^> 

one finds a set of equations of which the first two orders are 

5= - (1 + w)(6 - 3(f>) - 3H(c 2 s - w)S, (41) 

2 

9 = -H(l- 3w)9 — + -^k 2 S - k 2 a + k 2 cj), (42) 

1 + w 1 + w 

where c 2 = 5P/8p, and H = a/a. These equations can be recognised as the continuity and Euler equations 
respectively. 

Provided that the anisotropic stress term, a is negligible the hierarchy can be truncated after the 9 
term. While this is true for a perfect fluid it is not so for massive neutrinos. However, as discussed in 
detail in Rcf. 56 the anisotropic stress for massive neutrinos in the non-relativistic regime can to a good 
approximation be simply related to the c 2 term. The end result is a set of fluid equations for the massive 
neutrinos which is approximately 

6 = - 9, (43) 

'3 

—j 
2 



H9 - [\h 2 - %k 2 c 2 d ) 6, (44) 



where 



/ d 3 pf v 2 



-«t- f - . • ( 45 ) 

J d 6 pf 

For non-relativistic neutrinos the anisotropic stress term acts like an additional pressure term which 
prevents the growth of structure. Below a scale corresponding to 

Q H 2 
k 2 -—" (46) 

10 c d 
structures cannot form. For much smaller fc the solution corresponds to the growing mode of the CDM 
perturbations. In Fig. [2] I have plotted the free streaming scale as a function of a, as well as the comoving 
Hubble scale, aH. This clearly identifies the relatively narrow region in k where modes can grow. When 
a increases the neutrino velocity dispersion decreases and the region where growth is allowed becomes 
correspondingly larger. For a smaller neutrino mass the growing region would be smaller because fcps starts 
to increase with a only once neutrinos become non-relativistic. It is now also clear why the onset of growth 
of \I/ seen in Fig.[T]happens later for larger k. The given mode crosses into the region where growth is allowed 
at a later stage in cosmic evolution. 

We are thus able to identify three different regions for r and k: 

1) r < t(T = m): Neutrinos are relativistic and there is no growth of structure, 

2) t > t(T = to), k < fcps: structures grow as for CDM, 

3) t > t(T = to), k > fcps: Structures cannot form. 

When measuring fluctuations it is customary to use the power spectrum, P{k, r), defined as 

P(k,r) = \S\ 2 (r). (47) 

The power spectrum can be decomposed into a primordial part, Po(k), and a transfer function T 2 (k, r), 

P(k,T)=P (k)T 2 (k,r). (48) 

The transfer function at a particular time is found by solving the Boltzmann equation for S(t). 

In a universe with a mixture of neutrino hot dark matter and cold dark matter, the suppression of 
fluctuations on small scales can be found using an analytic approximation. On scales which are smaller than 
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Figure 2. The upper black line is the neutrino free streaming scale, fcpsi f° r a model with 5Zm„ = 1.2 eV. 
The lower black line is the comoving Hubble scale, aH. All modes in the shaded region can grow, whereas 
modes above the shaded region are subject to free streaming suppression. Figure adapted from 56 . 



the free-streaming scale for neutrinos, the usual Meszaros equation for the evolution of CDM perturbations 
is given by 



la--L(l-/„)J = 0, 

T T z 



(49) 



i.e. the gravitational source term is modified because neutrinos do not contribute to the gravitational 
potential, only to the background evolution. The normal growing solution 6 oc r 2 is then modified to 
5 oc r^v i+24(i-/„)-i)/2^ anc [ m j-jjg ijmit of /„ <C 1 the power spectrum is modified roughly by [I] 

AP 

-p 8/„. (50) 

The same suppression factor is found from an exact numerical solution to the Boltzmann hierarchy [78] . In 
Fig. |21 1 have plotted matter power spectra for various different neutrino masses in a flat ACDM universe 
(fl m + fl v + fl\ = 1). The parameters used were fib = 0.04, Hcdm = 0.26 — fl v , Ha — 0.7, h = 0.7, and 
n = 1 and the calculations were done with the CAMB [55 Boltzmann solver. 

The effect of massive neutrinos on structure formation only applies to the scales below the free-streaming 
length. For neutrinos with masses of several eV the free-streaming scale is smaller than the scales which can 
be probed using present CMB data and therefore the power spectrum suppression can be seen only in large 
scale structure data. On the other hand, neutrinos of sub-eV mass behave almost like a relativistic neutrino 
species for CMB considerations. The main effect of a small neutrino mass on the CMB is that it leads to 
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Figure 3. The matter power spectrum divided by the ACDM matter power spectrum for various different 
neutrino masses. The lines are for 5Zm„ = 0, 0.15, 0.6, and 1.2 eV in descending order. 

an enhanced early ISW effect. The reason is that the ratio of radiation to matter at recombination becomes 
larger because a sub-eV neutrino is still relativistic or semi-relativistic at recombination. 



4- 3. Non-linear evolution 

As was seen in the previous section the effect of neutrinos on the matter power spectrum in the linear regime 
can be reasonably well approximated by 



p — 



Pm{Qm = flc + Qu) for k -C fcFS 

(i - 8f u )p m {n m = n c + n u ) forfc>fc FS 



(51) 



with a smooth transition stretching over roughly 2 decades in k. 

However, at z = non- linear corrections become important already at k ~ 0.1 h/Mjpc. These non-linear 
corrections can be described either semi-analytically using the halo-model formalism, or using an extension 
of perturbation theory. Analytic calculations of the matter power spectrum with massive neutrinos included 
have been carried out in a number of recent studies [57J EH1 HH HO] • 

Ultimately, however, the validity of these approximations must be tested against N-body simulations. 
However, massive neutrinos pose a serious problem for N-body simulations since they have very high thermal 
velocities. For realistic masses most neutrinos will have much higher thermal velocities than the average 
gravitational streaming velocities, meaning that most neutrinos do not cluster in bound halos. Simply 
treating neutrinos as particles is difficult because the thermal motion introduces noise which completely 
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dominates the behaviour of neutrinos on small scales. However, it has been shown that the matter power 
spectrum can be computed at the 1% level of precision with neutrinos included, using N-body techniques 
[641 1651 l66l 167) . These simulations have shown a very interesting connection between thermal neutrino motion 
and gravitational streaming motion. In general neutrinos cause a suppression of matter fluctuations which is 
larger than the linear theory result, particularly on scales corresponding to the free-streaming scale at a given 
redshift, i.e. on scales where the virial velocity of a halo is comparable to the thermal velocity of a typical 
neutrino. The maximal suppression is roughly — 9.6/,,, peaked at scales around k ~ 0.5— 1 ft,/Mpc. Provided 
that other aspects of the power spectrum can be modeled sufficiently precisely this feature is a smoking gun 
signature for the presence of massive neutrinos. The presence of this neutrino induced depression in the 
power spectrum cannot be mimicked by other types of dark matter or dark energy, and furthermore it is 
sufficiently far removed from the baryon acoustic peak (at k ~ 0.15/i/Mpc) to be uniquely identified. It 
should also be noted that a peak suppression of ~ — 9.6/,, was also found in a recent study by Viel and 
Springel [63 which focussed mainly on simulations of the Lyman-a forest. Finally it should be noted that 
calculations using various extensions of linear perturbation theory, as described in [571 [SHI I5S1 ISO] show good 
agreement with simulations at high z and for relatively small k, i.e. the regime where nondinear corrections 
are mild, but not for larger values of k, i.e. k > 0.3/i/Mpc. 
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Figure 4. The matter power spectrum divided by the ACDM matter power spectrum for various different 
neutrino masses. The lines are for ~^2,ra v = 0, 0.15, 0.6, and 1.2 eV in descending order. The green lines are 
the linear theory predictions (figure reproduced from I119| V 



4-4- Other aspects of non-linear structure formation 

The matter power spectrum is not the only interesting large scale structure observable. Other probes which 
will be of importance in the next decade are for example weak gravitational lensing and the measurement of 
galaxy clusters. In [67] properties of dark matter halos in cosmologies with massive neutrinos were studied 
and found to be significantly different from standard ACDM halos. In general halos of a given mass form 
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later and are less concentrated in such models, and the number of halos of a given mass is also changed. 
Part of the effect comes from the different linear theory initial condition, but there is also a substantial 
effect from non-linear evolution. In |74j possible neutrino mass constraints from a future cluster survey were 
studied in a Fisher matrix analysis and using the Jenkins et al. [120] prescription for propagating the linear 
theory transfer function to the halo mass function. It was found that some of the cluster surveys which will 
be carried out in the next decade will allow for 0.05 eV sensitivity to neutrino mass. 

4--4-1. The halo mass function In |67j detailed studies were performed on how clusters form in cosmologies 
with massive neutrinos. The numerical set-up was based on the code described in |66j and from the N-body 
simulations, halos were identified using the AMIGA halo finder [75] . 

From the list of halos the number density of halos per mass interval, i.e. the halo mass function (HMF) 
can be constructed. In Fig. [5] the HMFs for cosmologies with different neutrino masses are shown. As 
expected the HMF is more suppressed in cosmologies with a larger neutrino mass. The suppression is largest 
for the heaviest, late forming halos. 

In Fig. the HMFs calculated from ./V-body simulations are compared with HMFs from the Sheth- 
Tormen (ST) semi-analytic formulae 80 . The ST fit is based on the fact that, as first pointed out by Press 
and Schechter [ST], the HMF can be written as 

MdM dn(M,z) „, .dv . . 

— diw-^ /M - (52) 

with v = [5 sc {z)/a(M)] 2 , where 8 sc (z) — 1.686 is the overdensity required for spherical collapse at z, and 
p = il m Pc- dn(M, z) is the number density of halos in the mass interval M to M + dM. The variance of the 
linear theory density field, a 2 (M), is given by 

2 f dk k 3 P Un (k) , 



a\M) = J - ^' \W{kR)\\ (53) 

where Pn n (k) is the linear theory matter power spectrum, and the Top-Hat window function is given by 
W(x) = (3/a; 2 )(sina; - xcosx) with R = {ZM/Airp) 1 ^ . 
The ST fit to vf{v) is 

"/M^(i + ^)(^)" 2 ^, (M) 

with v' — 0.707f and p = 0.3. A = 0.3222 is determined from the integral constraint J f^dv = 1. 

The upper panel in Fig. [6] shows that the agreement is poor if Q m = fi c + fli, + Q v is used in the ST 
formalism. However, this is due to a wrong definition of the halo mass: Even for the very largest cluster halos 
the neutrino component contributes very little to the halo mass. In reality, the mass inside the collapsing 
region should be calculated using J7 c + J7b, not Q, m . This amounts to neglecting the weakly clustering neutrino 
component when calculating the halo mass. The two lower panels in Fig. [B] shows the same ST fit, but using 
fl c + fib instead of J7 m . In this case the ST HMFs provide an excellent fit to the relative change to the HMF 
caused by neutrinos. As the figure at the bottom clearly demonstrates, the agreement is better than ~ 3% 
at halo mass scales where the ./V-body HMFs are accurate. Although the absolute HMFs, even for CDM 
simulations, do not match the ST HMFs more precisely than at the ~ 10% level, the relative change from 
adding neutrinos can be calculated significantly more accurately. 

Jj..Jf..2. Individual neutrino halos Detailed N-body simulations can also be used to probe the density profiles 
of individual halos. In [57] the overdensity of neutrinos in halos of varying size was studied for a variety 
of different neutrino masses. The results of these simulations are shown in Fig. [7J The figure shows three 
different types of curves: One for all halos of a given total mass, one with only isolated halos, i.e. halos which 
are not subhalos of a larger halo, and one which is calculated using the N-l-body method [73] (see [57] for 
details). 

The gravitational effect of a host halo is relatively much more important for neutrinos than for the CDM 
component: Due to free-streaming neutrinos will almost completely stream out of small halos (~ 10 12 M Q ), 
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Figure 5. Absolute and relative halo mass functions for 5 different neutrino cosmologies. The halo mass 
functions have been splined and smoothed together to obtain sufficient accuracy in the halo mass range 10 12 
to 10 15 Mq (figure reproduced from 67 ). 



and any measured value 8 U > will be caused by the host halo. The radial profile of b~ v will therefore be a 
superposition of a dominant flat profile from the host halo on top of a sub-dominant contribution from the 
~ 10 12 M Q halo itself. This fact can be seen in Fig. [7] 

In the same study it was also shown that the Tremaine-Gunn bound for neutrinos is in practise never 
saturated because neutrinos only start clustering very late in the evolution of the universe. In practise even 
the central density of neutrinos is significantly lower than what is estimated from the TG bound for a given 
neutrino mass. Finally, in Fig. [5] I show an example of a neutrino halo in a heavy host halo of 5 • 1O 14 M0. 
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Figure 6. Relative halo mass functions for different neutrino cosmologies compared with the predictions 
from the Sheth-Tormen formulae (black lines). Top: With Q m = Q c + SI;, + Q v in the ST formulae. Middle: 
With Q c + Qf, used instead of Q m in the ST formulae. Bottom: Differences between the TV-body and the 
ST predictions (figure reproduced from 67 ). 
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Figure 7. Neutrino halo profiles for Y2 m v = 0.3 eV (top), ~^2ra v = 0.6 eV (middle) and ^2,m u = 1.2 eV 
(bottom) for halo masses of 10 12 , 10 13 , 10 14 and 10 15 Mq. Profiles are calculated with the ./V-one-body 
method (dotted) and the Af-body method with a halo isolation criterion (solid) and without (dot-dashed) 
(figure reproduced from [57]). 
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Figure 8. CDM and ^2 rn v =1.2 eV neutrino distributions for a halo mass of ~ 5 ■ 10 14 Mq . Dimensions in 
each image is 5 h~ 1 Mpc. The images correspond to CDM, total neutrino, and q/T = 1 to 6 from top-left 
to bottom-right. Individual neutrino TV-body particles can be identified (figure reproduced from 67). 

5. Current constraints on the neutrino mass 

5.1. Parameter estimation methodology 

Because massive neutrinos affect structure formation it is possible to constrain their mass using a combination 
of cosmological data. The standard approach to cosmological parameter estimation is to use Bayesian 
statistics which provides a very simple method for incorporating prior information on parameters from other 
sources. Using the prior probability distribution it is then possible to calculate the posterior distribution and 
from that to derive confidence limits on parameters. There are standard packages such as CosmoMC [85], a 
likelihood calculator based on the Markov Chain Monte Carlo method [531 EH], available for this purpose. 

Although it is standard practise in cosmology to analyze data using a Bayesian framework there are 
examples of frequentist analyses. In most cases there is little difference in the inferred parameter ranges as 
long as the given parameters are well constrained by the data. 

Independent of the statistical method used results will in general depend on the number of parameters 
used to fit the data. Because of parameter degeneracies bounds on a given parameter will in general get 
weaker if more parameters are used in the fit. However, the question remains as to how many parameters 
should plausibly be included. Particularly with early data it was always a problem that changes in parameters 
like the neutrino mass could be mimicked by a combination of changes in other parameters, i.e. there were 
severe parameter degeneracies in the data. In such cases the actual constraint on a given parameter can 
depend strongly on the chosen model space. For example the constraint on neutrino mass from CMB and 
large scale structure data would degrade by a factor of two or more when the equation of state of dark energy 
was allowed to be different from -1. 

With current data this is less of an issue simply because better data causes many parameter degeneracies 
to be broken. In some cases this is true for just a single set of data, but in other cases the breaking of 
degeneracies arises from the combination of two or more data sets. 



5.2. Current bounds 

Unlike the first WMAP data, the WMAP-7 data release on its own provides a quite stringent constraint on 
the sum of neutrino masses of J2 m v < 1-3 eV at 95% c.l. [1 (down from ~ 2.1 eV in the first releases |84) ) . 
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However, this remains true only within the framework of the standard ACDM + m„ model. For example, 
even if a prior on Hq corresponding to the latest HST result [85] is added to the CMB data the bound 
worsens to Yl m ^ < 2-1 eV at 95% c.l. if N v ^ 3.04 is allowed and further to 2.6 if the dark energy equation 
of state, w, is allowed to be different from -1 [86 . 

Including large scale structure data the bound can be improved significantly. At present, by far the 
best large scale structure data comes from the Sloan Digital Sky Survey, for example in the form of the 
halo power spectrum (HPS) from the SDSS-LRG DR7 sample presented in [57], or from the BAO feature 
determined from the full SDSS-DR7 data set [88] . 

When combined with CMB data the HPS provides an upper bound of 0.61 eV in the minimal model 
(going down to 0.44 eV if the HST prior is imposed), and 1.16 eV in the extended ACDM+m„+N u +w 
model. The corresponding CMB plus BAO constraint is 0.85 eV for the minimal model and 1.4 eV for the 
extended model [55]. Note also the recent analysis performed in [5T] of the photometric SDSS MegaZ LRG 
DR7 catalogue which gives a very stringent bound of 0.68 eV when combined with other data. Ref. [62] 
contains an interesting analysis of the robustness of the mass bound against astrophysical uncertainties. 

Another similar SDSS related probe is the SDSS MaxBCG catalog, used in [55]. Effectively this provides 
a bound on the small scale normalisation of fluctuation power, and yields an upper bound on 0.4 eV for the 
minimal model. Another probe which in principle is very sensitive to the small scale power is the Lyman-a 
forest. The most recent and sophisticated analysis of Lyman-a data [63] gives an upper bound of 0.9 eV for 
the sum of neutrino masses. 

In conclusion the current bound on the sum of neutrino masses can be in the range between 0.3 and 
more than 2 eV, depending on the data and parameters used. In general: 

a) The bound on ^m^ for the minimal ACDM plus neutrino mass (8 parameters in total) is in the 0.4 eV 
range if CMB and LSS data is used. 

b) This bound can be relaxed somewhat when more parameters, such as w and N v , are included. In the 
most conservative case the bound is above 2.5 eV if only CMB data is used. 

c) When CMB data is combined with LSS data in the linear or almost linear regime, combined with a prior 
on the Hubble parameter the upper bound is robustly below 1 eV. This is true even for extended models. 

Here it should perhaps also be noted that the bound on neutrino mass from cosmic structure formation 
applies to any other, hypothetical particle species which decouples while still relativistic. This could for 
example be low mass sterile neutrinos, as discussed for example in [751 176] . It could also be relatively high 
mass axions which decouple after the QCD phase transition [55] |5T]. While there are some quantitative 
differences, the qualitative argument is the same for all such cases. In Fig. [9] I show one such example of a 
combined analysis of neutrinos and axions. 

6. Is there evidence for extra neutrino density? 

The WMAP data have consistently pointed to a value of N u larger than 3 as the best fit, and this is 
particularly true for the latest WMAP-7 data release [1 . With the WMAP-7 data plus the SDSS-LRG DR7 
data it was found by the WMAP team that N v = 4.25J^'gQ PQ, and similar results have been obtained in 
other similar analyses. While not statistically significant at more than 2<r it is nevertheless suggestive, and 
also leads to the question of what this extra energy density could be. 

In terms of neutrinos the simplest possibility by far is one or more light sterile states mixing with 
the active flavours. Such neutrinos can be either partly or completely thermalised via mixing in the early 
universe, with their final abundance depending strongly on the mass difference and mixing angle with the 
active species. 

For a fully or almost thermalised sterile state the usual cosmological mass bound applies and means 
that it cannot have a mass above 0.5-1 eV in order not to be in conflict with structure formation constraints 
[751 176] . However, for small mixing angles full thermalisation does not happen and the mass bound is 
significantly relaxed (76) . This would for example be the case for keV sterile neutrino warm dark matter. 
The most recent study of light sterile neutrinos in cosmology is the one performed in [77] . 
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Figure 9. 2D marginal 68% and 95% contours in the Y2m v -m a plane. The blue lines correspond to our 
old results using WMAP-5, the magenta lines our new results using CMB+HPS, and the green lines using 
CMB+HPS+HST (figure reproduced from |5T) ). 



Another possibility which is perhaps slightly less natural from a model building perspective, but not 
excluded, is the presence of large neutrino chemical potentials. In terms of energy density and non-zero 
chemical potential automatically leads to an excess energy density, with the relation between N v and £ = /i/T 
being 



N v = 3- 



£?' e 



7 



(55) 



taking into account the possibility of different chemical potentials for the thee flavours. 

Generically, large neutrino chemical potentials are excluded by the nucleosynthesis bound discussed 
previously. However, as was shown in 70 , by tuning the initial asymmetries carefully it is possible to 
circumvent the nucleosynthesis bound and achieve a significant contribution to N v from non-zero chemical 
potentials. 

Data from the Planck satellite which is currently observing the cosmic microwave background will have 
a much improved sensitivity to N v , probably reaching a{N v ) ~ 0.2 [}52]. If the best-fit value remains at 
N v ~ 4 Planck should therefore be able to provide evidence for new physics in this channel at fairly high 
significance. 
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7. Future probes of neutrino mass in cosmology 

In the coming years a variety of different new experiments will probe cosmic structure with high precision 
across a wide range in both scale and redshift. Some probes are based on refinements of techniques already 
used. Examples of this are the Planck mission which is currently measuring the CMB anisotropy and 
upcoming galaxy surveys such as BOSS. Others will use completely new techniques, for example looking at 
structures at high redshift using radio observations, or probe weak gravitational lensing of galaxies by cosmic 
large scale structure. 

In this bewildering multitude of different experiments there are some reasonably simple guidelines which 
can be used to estimate the sensitivity to neutrino parameters. The underlying quantity measured by most 
of these experiments is the matter power spectrum. When probing for example inflationary models higher 
order statistics such as the bispectrum or trispectrum can also be important, but for neutrino masses the 
power spectrum is normally the relevant quantity. 

The precision with which the power spectrum can be measured depends on a number of parameters. 
First of all, there is sample variance, i.e. the fact that only a finite number of Fourier modes are contained 
within the survey volume. In some cases (like the CMB) this is a fundamental problem related to the horizon, 
but in most cases it is a technical problem because the effective survey volume is determined by the size of 
the telescope, the time available, and the luminosity of the objects studied. Beating down sample variance 
is best done using bright source objects to trace the underlying matter distribution. However, this leads to 
another fundamental problem, shot noise. The density of bright galaxies is low and therefore small scale 
modes cannot be properly sampled. For galaxy surveys this is not a major concern because luminous red 
galaxies, such as used in the SDSS, are common enough that the power spectrum can be be measured to k 
significantly larger than 0.1 /i/Mpc. For weak lensing surveys it is more of a concern because many more 
galaxies need to be measured to infer an average distortion in a given patch of the sky. 

Finally, there can be serious systematic problems related to the extraction of the underlying matter 
power spectrum from a given measurement. In most cases the measurement is not of the matter density 
itself, but rather of some possibly biased tracer of the matter density. This is true for example in galaxy 
surveys where it is necessary to understand how the galaxy clustering power spectra for various types of 
galaxies are related to the underlying matter power spectrum. For Lyman-a forest measurements the problem 
is even worse because the measurement is of the column density of neutral hydrogen along lines of sight, 
a quantity which requires extensive modeling before a matter power spectrum estimate can be extracted. 
Similar problems will certainly face future measurements of temperature fluctuations in the 21-cm band. 

The following subsections contain a brief review of the techniques which are currently used or under 
development. Note that all sensitivities are to the sum of neutrino masses, ^m„, and at la unless stated 
otherwise. 

7.1. The cosmic microwave background 

While the WMAP satellite has provided cosmic variance limited measurements of the temperature anisotropy 
up to I ~ 500, there are vast improvements to be made. The Planck mission [92 is currently measuring the 
CMB anisotropy and will provide cosmic variance limited temperature measurements up to I ~ 2000, as well 
as vastly better measurements of the polarisation anisotropy. Even though the primary CMB signal is not 
highly sensitive to the neutrino mass, Planck will improve the WMAP constraint from 1.3 eV in the minimal 
model to 0.4-0.5 eV [H |53J |55], and to 0.5-0.6 eV in more extended models. Even though this result is on 
par with the current constraint it is significantly more robust because it relies on only one data set, which 
is furthermore measured on scales large enough to make non-linear contributions subdominant. 

However, even if such non-linear corrections are subdominant they are potentially important for 
measuring the neutrino mass. Just as with any other source of radiation at high redshift, the CMB signal 
is affected by gravitational lensing along the line of sight. For the CMB this effect can be diagnosed in a 
number of ways. For example it introduces non-Gaussian features in the otherwise Gaussian primary signal 
(see e.g. [9" H l93]). Another important effect is that it mixes E and B mode polarisation of the CMB [9"5"ll96| . 
and any experiment with sufficient sensitivity to the very weak B mode signal should be able to identify 
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this. 

With lensing extraction Planck data alone should be able to reach a sensitivity of 0.1-0.15 eV |93j . 
depending on the complexity of the model assumed. Further into the future, experiments highly sensitive to 
B-mode polarisation such as CMBPol [97] might plausible reach a sensitivity of 0.04 eV. 

7.2. Galaxy surveys 

Galaxy surveys have up to now been the most direct way of measuring the matter power spectrum on 
intermediate and small scales, and therefore also the most direct probe of the suppression of fluctuation 
power caused by the presence of massive neutrinos. At present by far the largest spectroscopic survey is the 
SDSS, and, together with the WMAP CMB data, it provides an upper bound of approximately 0.6 eV on 
Em, [57]. 

Galaxy redshift surveys measure the power spectrum of galaxy number density fluctuations, P g {k). In 
turn, this power spectrum is related to the underlying matter power spectrum P{k) via 

P B (k) = b 2 (k)P(k), (56) 

where the bias parameter b depends on both scale and on the type of galaxies surveyed. This has been shown 
to be a significant problem on some surveys and therefore emphasis has shifted towards basing surveys on 
luminous red cluster galaxies which constitute a fairly homogeneous sample. This is for example the case 
for the SDSS-LRG sample which was used to derive the current 0.6 eV upper bound. 

A number of larger galaxy surveys will be carried out within the next decade and will increase sensitivity 
to neutrino mass significantly. Given some survey design one can expect to measure the galaxy power 
spectrum P g {k) up to a statistical uncertainty of 98 



AP g (k)= '" 



2tt w(k) Alnfc 



tin 



(57) 



Here, w(k) — (fc/27r) 3 V c g , V e g is the effective volume of the survey, n g is the galaxy number density, and 
Aln/j is the bin size at k in lnfc-space. Future surveys will go deeper (and therefore have larger n g ) and 
wider (and therefore increased V c s) leading to much smaller errors on the power spectrum. 

The precision with which the power spectrum can, in principle, be measured is related to the survey 
volume because that is a measure of the number of independent Fourier modes available. On small scales 
precision is limited by shot noise, i.e. by the sparsity of galaxies. However, in practice this is not the most 
significant problem on small scales. Rather, the usefulness of small scale data is limited by the fact that 
structures are non- linear. At z = this effectively cuts away all data at k > 0.1 h/Mpc. 

However, most of the upcoming surveys aim at measuring at higher redshift than the SDSS and therefore 
the problem of non-linearity will be somewhat alleviated. In (99) a study of neutrino mass constraints was 



carried out for a number of proposed surveys combined with Planck data. Very roughly, the HETDEX [100] , 
WFMOS [TUB] , or BOSS [TUT] surveys, together with Planck should push the sensitivity to about 0.2 eV at 
95% C.L., and a future space-based mission such as JDEM or EUCLID could yield a sensitivity of around 
0.1 eV (95% C.L.). 

The major theoretical hurdles that need to be addressed in order to extract these sensitive limits are 
understanding the nonlinearities and bias. Simulations and cross-correlating with lensing surveys can help 
with these issues. 

7.3. Weak lensing surveys 

Weak gravitational lensing provides the only direct way to probe the cosmic matter distribution, as opposed 
to the distribution of some possibly biased tracer like galaxies. Perturbations in the matter density bend 
light rays from distant sources and leads to distortions in the shapes of objects such as galaxies. While 
galaxies are intrinsically elliptical the alignment of their axes should be close to random. Lensing on the 
other hand produced a coherent effect and by measuring the shape of many source galaxies and subsequently 
averaging it is possible to get an estimate of the weak lensing distortion field. 
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By measuring the weak lensing distortion (angular) power spectrum it is possible to probe the power 
spectrum of matter fluctuations between source and observer. Since it is at least in principle possible to 
measure the redshift of source galaxies it is furthermore possible to perform weak lensing tomography and 
probe the time evolution of the cosmic density field. It is impossible to obtain spectral redshifts for all 
source galaxies (the largest surveys will contain billions of such galaxies), but redshifts can be measured 
quite accurately by observing in different wavelength bands. 

Although weak lensing surveys are still in their infancy they already provide some constraints on 
cosmology. At this stage the constraint is mainly just on the amplitude of matter fluctuations on intermediate 
scales. Data from the CFHT Legacy Survey has been used in combination with the WMAP 5-year data to 
provide an upper bound on the neutrino mass of Yl m v < 1-1 at 95% C.L. J102J . 

Future weak lensing surveys from instruments such as the LSST [103] will cover a large fraction of the 
sky and extend to redshifts of 1 or beyond. When combined with Planck data this will allow for very tight 
constraints on the neutrino mass. For example |107j (see also |108j ) has estimated that by binning the LSST 
data in 5 tomographic bins a sensitivity of 0.07 eV at 95% C.L. can be reached. Proposed future surveys 
from space by experiments such as DUNE/Euclid may reach a comparable or better sensitivity [104] . 

One possibly severe problem with weak lensing surveys reaching this sensitivity is that, like CMB 
measurements, they measure the correlation function on a sphere, not in 3D. This means that they are 
sensitive to an integral along the line of sight of the underlying 3D fluctuation spectrum. This integration 
mixes small and large scale modes and therefore non-linear effects become important for weak lensing surveys 
even at relatively modest /-values. Since cosmological constraints can mainly be derived from moderate and 
large Z-values it means that non-linear corrections must be under control at the percent level. For neutrinos 
this has already been achieved [64 , 65 , 66 , but for example baryonic physics can contribute an uncertainty 
of 10% or more on scale which are relevant for cosmological parameter extraction |105j . 

7.4- Lyman-a forest measurements 

An indirect way of measuring the matter power spectrum at high redshift consists of measuring absorption 
lines from neutral hydrogen along the lines of sight to quasars. The width and amplitude of the Lyman-a 
absorption lines are related to the local hydrogen density and therefore also (albeit indirectly) related to 
the local matter density. An advantage of this method is that it can be used to probe high redshifts where 
structure are more linear. It can also probe much smaller scales than conventional galaxy surveys because 
it does not rely on a high density of bright objects such as LRGs. Historically the method was first used to 
constrain neutrino masses in |109j , and subsequently the large SDSS sample of absorption spectra was used 
in [HOj to infer an extremely stringent bound on J2 m » °f 0.17 e ^ a * 95% C.L. 

However, as has been discussed at length in |1101 II 1 11 155] the bound comes mainly from the fact that 



the small scale matter power spectrum normalisation is inferred in (110) to be significantly higher than the 
one found from other cosmological data. Using more up to date simulations which include neutrinos in 
the simulations [63] find a normalisation consistent with other probes and therefore a much less restrictive 
bound. 

In the future there will be a number of new surveys such as BOSS [10 lj which will probe the Lyman-a 
forest in more detail. Together with more accurate modeling it seems likely that such surveys can reach a 
sensitivity in the 0.1 eV region when combined with CMB data. 

7.5. Cluster surveys 

Galaxy clusters, with masses around 10 14 — 1O 15 M are the largest gravitationally bound objects in the 
universe. Any change in the matter power spectrum on small scales shifts the mass of the largest object 
capable of forming before the onset of dark energy domination, and since neutrinos have a big impact on the 
power spectrum on these scales the presence of massive neutrinos has a noticeable impact on the halo mass 
function, i.e. the number density of halos with a given mass. This point was discussed in detail in Section 
14.4.11 (and for even more detail the reader is referred to [67] ). 
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The biggest problem facing cluster surveys is that it is difficult to estimate the mass of a cluster based 
on simple observables such as their X-ray luminosity. Weak lensing measurements are more robust, but 
also significantly more time consuming. Fortunately, a number of large scale cluster surveys will become 
available in the not too distant future. For example the LSST survey [112] will provide a sample of thousands 
of clusters. Hopefully this will also make it possible to make a much better calibration of such quantities as 
the mass-temperature relation. 

If these issues can be reliably addressed, cluster surveys in combination with CMB data from Planck 
could potentially reach a sensitivity in the 0.04-0.07 eV range [74] . 

7.6. 21-cm measurements 

Radio band observations of the 21-cm line of neutral hydrogen can potentially reach much higher redshifts 
than any observations in the visible or infrared bands, and may provide access to information about the 
universe prior to reionization at 6 < z < 12. The brightness of a given patch of the universe in the 21-cm 
band depends on a number of factors and is given approximately by [113j 

where xh is the fraction of neutral hydrogen, S is the density inhomogeneity, T s is the spin temperature, 
and v r is the radial velocity. Thus, a measurement of 5T%\ in principle provides information on 5 and can 
be used to extract the matter power spectrum. 

One huge advantage of the 21-cm radiation is that it covers an incredibly large volume, comparable to 
the entire Hubble volume. Furthermore, since it is a 3D measurement, as opposed to the CMB measurements 
which are on a sphere, there are many more independent Fourier modes than in any CMB measurements. 

Given that the power spectrum can be reliably extracted from future observations, some of the upcoming 
or proposed radio experiments will be extremely sensitive to neutrino masses. For example, the Square 
Kilometer Array (SKA) might constrain ^2 tn v to 0.02 eV, while the proposed Fast Fourier Transform 
Telescope (FFTT) [111 could achieve a{J2m„) ~ 0.0003 eV [TTSlfTB] . 

However, before this can be achieved there are some severe obstacles to be dealt with. The foreground 



contamination of highly redshiftet 21-cm radiation is enormous |113j . For example the galactic synchrotron 
emission is 4-5 orders of magnitude larger than the primary signal. There are additional problems with 
terrestrial contamination. It remains to be seen to which extent (if indeed any) 21-cm measurements can be 
used for precision cosmology. 

7.6.1. Summary of future probes As can be gauged from the above discussion there is cause for optimism 
on cosmological neutrino mass measurements. Many of the surveys discussed above will be realised within 
the next 5-7 years, and a number of different techniques with completely different systematics all have 
comparable potential sensitivities in the o~(m v ) ~ 0.05 eV range. This means that there will be many cross- 
checks available and that any possible detection of a neutrino mass will not rely on any single data set which 
could be dominated by unaccounted for systematics. Table Q] shows a summary of the achievable sensitivities 
using various techniques. 

Another question is what exactly will be learned about neutrino physics from these measurements. If 
neutrinos hierarchical or inverted hierarchical the sum of neutrino masses will be in the 0.05-0.1 eV range, 
and cosmological probes within the next decade will at most be able to provide a l-2tr evidence for non- 
zero mass. This also means that it is unlikely that cosmology will determine whether the normal or the 
inverted hierarchy is the correct model in that time frame. In order for such a determination, the formal 
sensitivity would have to be significantly better than 0.05 eV and so far 21-cm observations and very large 
scale lensing surveys seem the best path forward (see |117l 1118] for a more detailed discussion of probing 
the hierarchy with cosmology) . If neutrino masses are hierarchical or inverted hierarchical it is also unlikely 
that direct laboratory measurements will provide any significant evidence of neutrino mass within the next 
decade. Certainly beta decay experiments like KATRIN will remain far off the mark for the foreseeable 
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future, whereas neutrinoless double beta decay experiments may or may not find evidence, depending on the 
nature of the hierarchy and exact Majorana phase structure. 

However, it is also entirely possible that neutrino masses are degenerate in which case cosmology should 
provide a high significance detection within the next decade. Furthermore this possibility is interesting 
because both beta decay and neutrinoless double beta decay experiments hold the potential to achieve a 
positive detection. Having measurements of three different new neutrino mass observables, m Ve , mpp, and 
y^ m„ will provide a plethora of new information on neutrino physics, including the values of Majorana 
phases. 



Probe 


Potential sensitivity (short term) 


Potential sensitivity (long term) 


CMB 


0.4-0.6 


0.4 


CMB with lensing 


0.1-0.15 


0.04 


CMB + Galaxy Distribution 


0.2 


0.05-0.1 


CMB + Lensing of Galaxies 


0.1 


0.03-0.04 


CMB + Lyman-a 


0.1-0.2 


Unknown 


CMB + Galaxy Clusters 


- 


0.05 


CMB + 21 cm 


- 


0.0003-0.1 



Table 1. Future probes of neutrino mass, as well as their projected sensitivity to neutrino mass. Sensitivity 
in the short term means achievable in approximately 5-7 years, while long term means 7-15 years. 



8. Detection of the cosmic neutrino background 

Standard model physics likewise predicts the presence of a Cosmic Neutrino Background (Ci^B) with a well 
defined temperature of T v ~ (4/ll) 1 / 3 T , 7 . While it remains undetected in direct experiments, the presence 
of the CvQ is strongly hinted at in CMB data. The homogeneous C^B component has been detected at 
the 4-5<7 level in the WMAP data. Furthermore, this component is known to be free-streaming, i.e. to 
have an anisotropic stress component consistent with what is expected from standard model neutrinos (see 
[mi H1II23IIM1 1123 H3 [113). Finally the standard model neutrino decoupling history is also confirmed 
by Big Bang Nucleosynthesis (BBN), the outcome of which depends on both the energy density and flavour 
composition of the C^B. 

While this indirect evidence for the presence of a Ci^B is important, a direct detection remains an 
intriguing, but almost impossible goal. The most credible proposed method is to look for a peak in beta 
decay spectra related to neutrino absorption from the C^B |128l 11291 1130] . although many other possibilities 
have been discussed [I5H H521 IT531 1731 H331 H551 WM H57) ■ 

Near the electron endpoint energy there should be a small contamination from anti-neutrino capture 
from the background. This process is monoenergetic, and the electron energy should be E e ~ Q + m v . 
Thus, the separation between the end of the beta spectrum and the absorption peak is 2m v . Depending on 
the energy resolution and the source luminosity this may or may not be visible in future experiments. The 
problem was recently studied in more detail in |129j , with the conclusion that currently planned experiments 
like KATRIN are still several orders of magnitude away from a detection. Nevertheless, the method is in 
principle feasible and so far seems the most promising way to detect the neutrino background. 

The neutrino absorption method was first investigated by Weinberg |128j , based on the possibility that 
the primordial neutrino density could be orders of magnitude higher than normally assumed due to the 
presence of a large chemical potential. Although a large chemical potential has been ruled out because it is 
in conflict with BBN and CMB, the method may still work and recently there has been renewed interest in 
detecting the C^B using beta unstable nuclei. 

Although the direct detection of the Ci/B is already very challenging, one might speculate on the 
possibility that in the more distant future anisotropics in the C^B will be detectable. 
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Figure 10. Primary C^B spectrum for different neutrino masses. 



8.1. Anisotropics in the cosmic neutrino background 

In the more distant future it may become possible to not only detect the cosmic neutrino background, but also 
to measure the equivalent of the CMB anisotropy. Interestingly, for massive particles the cosmic background 
anisotropy becomes markedly different from the massless case. This was studied in |119| for the case of light 
neutrinos. For neutrinos with average thermal velocity larger than v ~ 1000 km/s, very few neutrinos are 
bound, even in massive halos. Therefore, for neutrino masses below 0.1 eV a linear theory calculation of the 
cosmic neutrino background is appropriate. Furthermore it reveals some interesting features of how radiation 
anisotropies are generated. Since neutrinos are non-relativistic during a significant epoch in the late-time 
universe the gravitational source terms are much more important than for photons. Since these terms enter 
only through the lowest multipoles the effect is to very strongly enhance the low-/ part of the spectrum, an 
effect which can be seen in Fig. 1101 If it ever becomes possible to measure the neutrino background, one 
would therefore mainly expect a low-/ signal associated with the most massive mass eigenstate. 



9. Discussion 

I have reviewed the topic of neutrino cosmology. The topic is quite diverse and I chose to put particular 
emphasis on how neutrinos affect structure formation in the universe while other interesting topics such 
as leptogenesis were left out. At present neutrino physics presents perhaps the best example of how 
precision cosmology can be used to address particle questions which would normally be probed in laboratory 
experiments. The main example is the absolute mass of neutrinos: Absolute neutrino masses can only be 
measured at great difficulty in the lab, using either beta decay or neutrinoless double beta decay experiments. 
On the other hand, even small neutrino masses have a great impact on structure formation, and a sum of 
neutrino masses in the eV range would have produced a measurable suppression of small scale structure. 
The current upper bound on the absolute neutrino mass is approximately ^2 m v ^0.5 eV, depending 
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somewhat on the cosmological model and the data used. However, the bound cannot be pushed much beyond 
1 eV even if several additional cosmological parameters are introduced. In the future, a large number of new 
surveys will bring down the sensitivity to neutrino mass by an order of magnitude or more. For example 
weak gravitational lensing surveys will reach a sensitivity below 0.05 eV within the next decade. Even this 
precision may not, however, be enough to guarantee a detection if neutrino masses are strongly hierarchical 
and follow the normal hierarchy. To go below 0.05 eV at high significance will likely require future 21-cm 
surveys or possibly a new generation of space based weak lensing measurements. The timescale for such 
measurements is likely to be 15-20 years or more. However, the theoretical lower limit to cosmological 
neutrino mass measurements is extremely low, and certainly low enough for a measurement of the individual 
neutrino masses even in the normal hierarchy. At this stage it seems likely that the first positive detection 
of a non-zero neutrino mass (other than from oscillation experiments) will come from cosmology, and in any 
case cosmology will remain a powerful laboratory for neutrino physics which is complementary to laboratory 
experiments. 
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